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The stochastic model of classical system of particles (partons), which dynamics 
includes random walk in plane as well as processes of death, splitting, annihilation 
and fusion of partons, is considered. A set of equations for multiparticle distribu¬ 
tion functions for this system can be described in terms of diagrams of the Reggeon 
Field Theory (RFT) with supercritical pomeron, where rapidity plays a role of the 
time variable. The multiparticle inclusive distributions of partons correspond to 
multi-pomeron states in this analogy. In order to calculate the amplitude of hadron 
interaction at given energy we define an operator of linkage of any pair of partons 
from two parton sets (the projectile and the target). The form of the hadron- 
interaction (parton-linkage) operator is determined from the requirement for the 
result to be in correspondence to the RFT formulas. It is shown that the require¬ 
ment of Lorentz invariance of the amplitude is fulfilled in this interpretation only 
at definite relation between the probabilities of parton splitting, annihilation and 
fusion. Interpretation of various methods of s-channel unitarization is discussed in 
terms of the model. 
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Yet say not sadly: they have left us! 

But say, with gratitude: they were.” 

— Vasily Zhukovsky 

This paper is dedicated to memory of Misha Marinov. One of the first pa¬ 
pers published by Misha was a work done together with Roginsky on relation 
between s- and t-channel helicity amplitudes which was important for descrip¬ 
tion of spin particle interactions in Regge theory. I’ve just come to this field 
and was strongly impressed with mathematical elegance of this paper, espe¬ 
cially with the fact that the direction to the unphysical point (centre-of-mass 
of t channel) turned out to be real and physical. Unfortunately, this paper 
contained a minor mistake (change of helicity sign for transition from particle 
to antiparticle was not taken into account). I think that this mistake was re¬ 
lated with Misha’s style because it could be easily detected by consideration 
of the simplest example of irN scattering. But Misha liked general considera¬ 
tions so much that didn’t carry on elementary test. Later on these results were 
independently reproduced (without the mistake) in the well known paper by 
Trueman and Wick where examples were carefully considered. 

In Theoretical Department Misha acted as a receptor providing an in¬ 
formation about new physical ideas especially using advanced mathematical 
technique. He was in earnest about this mission and prepared these talks very 
carefully. He always had a detailed plan which he followed very hard and didn’t 
allow uncontrolled discussions and digressions if they were not planned. Due 
to this hardness he always managed to convey to the audience main ideas and 
his own relation to them. Up to now my attitude to many of these topics is 
affected by primary Misha’s interpretation. 

Misha was born Teacher, and it was a mystery for me why he had no 
students in ITEP to teach or to supervise them. Once I asked him to give a 
couple of popular lectures on path integral for me and some colleagues. He 
agreed instantly, and as a result we got a half-year course of lectures, which 
was transformed eventually to the remarkable review in Physics Reports. And 
how many such opportunities we missed! 

Everything what he did, Misha tried to execute professionally, either it 
was scientific paper writing, tea making, or karate training. I remember as 
during one of ITEP Winter Schools we were playing charades trying to show 
meaning of word or its part by performing short sketches. Misha chose to show 
just the word ’’charade” (in russian). First his team performed three sketches 
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for its parts, and then, to show the word as a whole, they repeated all these 
sketches again at a rapid pace, representing the process of the game. It was 
well done, and we could see perfectionism of Misha as a theater director. To be 
more professional as an actor he used tooth-paste for his make-up (performing 
a role of devil). 

Misha’s unlimited erudition and knowledge were always at service of his 
colleagues. They consulted with him on any subject and always got thoughtful 
discussion and useful references, often to old and forgotten papers. The pa¬ 
per which I present in this volume just illustrates the lack of Misha’s advices 
and consultations for me - when it has been finished I met the old paper by 
Grassberger and Sundermeyer where the same model for reggeon field theory 
was formulated 

I liked to discuss with Misha Jewish history, origin of Jewish names, prob¬ 
lems of Jewish young people at entrance examinations in universities, etc. But 
not less amount of exciting information I got from him on Russian history and 
language, literature, biology, and many other things (by the way, the origin of 
the word ” eikonal” appearing in this paper and common with the word ” icon” 
I knew from Misha). 

Misha’s high standards exerted influence on attitude to different situations, 
including ethical problems, in the circle of ITEP theorists. His opinion often 
was crucial giving, as a tuning fork, a point of reference. And now, what 
unphysical point Misha would be at, the orientation to his standards is still 
fixed for us. 


1 Introduction 

Processes of strong interactions at high energies are described usually in frame¬ 
work of reggeon theory EJ. Systematic description of high-energy interactions 
of hadrons in terms of reggeon field theory (RET) was introduced by Gribov in 
ref.B. These results were obtained in a hybrid field theory model by calculation 
of Feynman diagrams containing the reggeon amplitudes as building blocks. 

Experimental data on hadron interactions show that the total cross sec¬ 
tion is dominated by multiparticle processes and the number of secondaries 
increases with energy growth. Therefore only models taking into account large 
number of effective degrees of freedom (referred below as to partons for brevity) 
in a fast hadron can correspond to the nature of reggeon amplitudes. 

Attempts to give the parton interpretation of the reggeon amplitudes 
started immediately after appearance of the notion of the reggeon. In ref. S 


I am grateful to A. Kaidalov and M. Ryskin who drawed my attention to this paper. 
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a representation of the pole reggeon diagram was introduced as a multipe¬ 
ripheral fluctuation of the fast hadron. Simultaneous interactions of several 
multiperipheral fluctuations with the target correspond to nonplariar reggeon 
diagrams and explain the presence of the regge-cut contributionsu Later on 
interpretations of reggeon interactions were suggested using as inner degrees 
of freedom the objects of quotum chromodynamics - the quark-gluon string 
model 0, the BFKL pomeronQ, etc. 

The common feature of all these models is consideration of a fast mov¬ 
ing hadron as a complicated composite system which components (partons) 
are continuously distributed in rapidity values and randomly in the impact 
parameter plane. Hence, the Fock wave function of the fast hadron is a super¬ 
position of multiparton states and its structure depends on the hadron energy. 

Interaction of the hadron with a target is realized only due to interactions 
of the most slow (in the target frame) partons produced step by step from faster 
components (short-range interaction in the rapidity space). The amplitude of 
interaction is determined by the number of slow partons n{y) in the vicinity 
of the target. The probability of interaction with a target of only one of slow 
partons is proportional to n{y) and corresponds to the regge-pole contribution. 
Interactions of two and more partons correspond to the regge cuts of second and 
higher orders. The cut contributions are sign-alternating because of shadowing 
effects. 

Growth of hadron cross sections with energy gives evidences in favor of 
the so called supercritical pomeron (Ap = ap(0) — I > 0). RFT with the 
supercritical pomeron encounters a number of difficulties which are not resolved 
in satisfactory way. The pole amplitude corresponding to the one-pomeron 
exchange increases with energy as s^ violating s-channel unitarity. Diagrams 
of higher orders, in principle, have to restore the unitarity. However, when the 
energy increases the number of essential diagrams grows drastically, and there 
is no satisfactory way to sum a series of reggeon diagrams at present. Usually 
one takes into account only a set of the simplest diagrams corresponding to 
the eikonal or “quasi-eikonal” approximation which corresponds to account for 
the non-enhanced cut only and disregards inter-pomeron interactions. The 
scattering amplitude in this approximation looks as for scattering on almost 
black disc of radius proportional to the rapidity y in correspondence to the 
Froissart regime. In some approximations contributions of±he enhanced regge 
cuts were summed. E.g. in the Schwimmer approximationLI, which is justified 
for hadron-nucleus interactions, the fan-type enhanced diagrams are taken into 
account but reggeon loops are r^dpjted. Analysis of more complicated sets of 
diagrams was carried out in refsSfl’t3. However, the problem of construction of 
hadron scattering amplitude satisfying both s-channel and t-channel unitarity 
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is not solved up to now. 

In the parton picture the increase of the one-pomeron contribution with 
energy in the regime of the supercritical pomeron means that the number of 
soft partons grows exponentially with hadron rapidity, ri{y) ~ exp(A 2 /). Thus, 
the mechanism of parton multiplication has to be present in the course of 
the rapidity evolution of the partonic fluctuation. It could be splitting of the 
multiperipheral ladder, gluon cascading, etc. From general arguments fusion 
of partons may also take place beside splitting. 

The kinetic model of parton interactions equivalent to RFT was suggested 
in paper by Grassberger and Sundermeyer EiJ. In fact, this model gives a 
desription of stationary structure of the Fock wave function of a fast hadron. 

In presented paper we discuss interaction of two parton systems in this 
model We analyze a possibility of consistent probabilistic interpretation 
of the interaction amplitude in this approach and, particularly, its Lorentz 
invariance. 

The main assumption of the model is continuous parton evolution in rapid¬ 
ity, i.e. an origin of slower partons from the faster predecessors. In the course 
of evolution partons undergo a random walk in the impact parameter plane. 
Other essential ingredients of parton dynamics are processes of splitting of one 
parton into two ones and fusion of two partons into one parton. For general¬ 
ity we consider also vanishing (“death”) one of partons at given rapidity and 
“annihilation” of pair of partons at close points. The correspondence can be 
proved between multiparticle parton distribution functions and multipomeron 
vertices of the RFT with lagrangian including the elementary three-pomeron 
coupling and the pomeron-scattering term. Thus, the pomeron in this model 
is not some specific object but just one-parton inclusive distribution, while 
many-pomeron states are connected to the inclusive multiparton distributions. 

In order to introduce into the model hadron-hadron interactions it is nec¬ 
essary to define the operator of the interaction of two parton systems with 
different rapidities. We show here that this operator can be defined in the way 
maintaining the correspondence with the interaction amplitude given by sum¬ 
ming all diagrams of RFT including loop diagrams. It is shown that if some 
definite relation between the constants of splitting, annihilation and fusion of 
partons holds, the hadron interaction amplitude depends only on the sum of 
rapidities of the interacting hadrons according to the Lorentz invariance re¬ 
quirement of the RFT. However the price for the correspondence with RFT 
is a loss of the probabilistic interpretation, because the interaction operator 
doesn’t have, in general, a probabilistic meaning. Nevertheless, the interaction 

^ Unliinkilv. T found thp nanerFl only after this text has been finished. I elaborated the same 
kinetic model for RFT so the content of first section overlap considerably with ref. EdI. 
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amplitude does not violate s-channel unitarity, if initial parton distributions 
can be represented as superpositions of Poisson distributions. 

The layout of the paper is as follows. In sect.2 the model is formulated 
in simplified version which does not take into account parton diffusion in the 
transverse plane. In sect.3 the diffusion of partons in the impact parameter 
plane is included and equations for the multiparticle distribution functions 
are derived both in coordinate and momentum representations. The equiv¬ 
alence of the equations for these multiparton distributions and equations for 
the many-pomeron vertices in the “cut” RFT is demonstrated in sect.4. In 
sect.5 we define the interaction operator for two parton sets and demonstrate 
independence of inelastic cross section on a choice of Lorentz frame. In sect.6 
the main model approximations in RFT are considered — the eikonal approx¬ 
imation, the Schwimmer model and some others. The picture of an expanding 
black disc corresponding to the Froissart regime is discussed. The problem 
of s-channel unitarity is discussed in sect.7. Finally, in sect.8, we summarize 
the main results of the paper and discuss briefly possible generalizations and 
applications of this approach. 


2 The model without diffusion 

We discuss in this section the 0-dimensional model (without parton diffusion) 
which posesses some typical features of the general case. This approximation 
becomes adequate when it is possible to ignore parton diffusion, e.g. in the 
regime of the uniform density distribution or in the case of interactions with 
heavy nuclei. 

Let us consider a system of partons with variable number of particles N, 
which evolves in time y due to the following elementary processes: death of 
a parton, annihilation of two partons, splitting of any parton into two ones 
and fusion of arbitrary parton pair into one parton. Denote the probability of 
the iV-parton state at the moment y as Pn(j/), and the probabilities of death, 
annihilation, splitting and fusion of partons for the unit of time as mi, m 2 , A 
and v, respectively. 

The evolution equations have the form: 


= -(A -I- mi)Nptf - {v + m 2 )N{N -l)pf^ + X{N - 1 )pn-i 
ay 

+ mi{N + 1)pjv+i + + l)iVpAr+i + m 2 {N + 2){N + 1)pn+2 , > 2) , 

dpo , „ 

— = mi Pi -I- 2m2 P2 , 
ay 



Probabilistic Model of Reggeon Field Theory 7 


—— =-{X + mi)pi+2{mi + v)p2 + 6m2P3 ■ ( 1 ) 

ay 

The initial condition for eqs (|^ is given by the parton distribution at the 
moment y = 0, and the conservation of the total probability p(y) = PNiv) 
follows from (|^): 



( 2 ) 


The equivalent set of equations can be written down for the factorial mo¬ 
ments ps, defined as expectations of quantities {N)s = N!/(N-s)!: 


Ps=^{N)sPn=^ , (s = 0,1,2,...) . (3) 

N=0 N=0 ^ 

After multiplying the both sides of the eqs ( 0 ) to {N)s and expressing the 
products appearing in the right-hand sides in terms of {N)s, (A^)s_i, ... , we 
come to the set of equations for the moments: 


dps 

dy 


As(s - l)ps-i -I- (A - mi)sps 
— (i/ + m2)s{s — l)ps — + 2m2)sps+i ■ 


(4) 


In particular, the equation for p{y) = po{y) is reproduced at s = 0. 

It is convenient to define the generating function G(w\y), which is at the 
same time the exponential generating function for moments: 


CX5 OO / 1 A 

G{w]y) = ^ PN{y)w^ = '^Ts{y )^^-^^— 

N=0 s=0 


(5) 


It satisfies the partial differential equation: 

r)G r)G rP G 

— = (1 - w){mi - Xw)- -h (1 - w)[m 2 + {m 2 + v)w\—-^ 

oy aw ow‘‘ 

with initial condition 


= Go{w) = ^ Pn{0)w^ . 

N^O 

The boundary condition at rc = 1 is maintained automatically: 


( 6 ) 

(7) 


G(l;y) = l . 


(8) 
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In case mi = m 2 = 0 it holds at w = 0: 

G{0-y) = Go{0) . (9) 

Let us consider some special cases: 

(a) The absence of parton fusion and annihilation (v = 0, m 2 = Oj 

In this case Eq.(||) becomes of 1st order and after the change of variables 
w ^ 7 , where 


d'y = 


dw 


Xw^ — (A + mi)w + mi 
one obtains 

dG{r.y) dG{r,y) 


, , ( 10 ) 


1 — exp[(A — mi) 7 ] 


dy dj 

G(7;0) = Go(7) , G(-oo;2/) = 1 


( 11 ) 

( 12 ) 


with notation G( 7 ; y) = G{w{'y)] y), Go( 7 ) = ^' 0 (^( 7 )). Its solution is a func¬ 
tion of the sum of the variables j + y satisfying the initial condition (Q), 


G( 7 ;y) = Go(7 + 2/) , 

or, coming back to the variable w, 

G{y, w)=J2 W{y)^pN{0) = Go{W{y)) 


(13) 

(14) 


N 


where 


Wiy) = 


1 — r] w — T] — ri{w — 1 ) 


1 — e^G+y) w — r] — {w — 1) e^y 
A = A — mi , rj = mi/A . 


(15) 


This distribution differs significantly from the Poisson distribution being char¬ 
acterized by a presence of strong correlations. 

If there is Nq particles in the initial state, then 


G{w,y) = [W{y)] 


No 


(16) 
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The mean value of particles increases exponentially with energy in this case: 

Ml( 2 /)=iVoe^^ (17) 


and the second moment behaves as follows: 

M 2 ( 2 /) = IVo [iVo + . (18) 

V 1-^7/ l-?7 


In case when toi ^ 0, the state without partons {N 
absorbing state. The probability po( 2 /) can only increase: 


Po{y) = 


'q{e^y — 1 ) 

gAy _ 


0 ) is so called 


(19) 


If the process of parton death dominates (A = A — toi < 0), then pa{y) 1 
at y ^ oo. If A > 0, this probability goes to the constant rj = mi/X. 

In the case of mi = 0 the probability po{y) = 0, and the probability 
distribution can be written as 


yiV —JMo j 


(b) The asymptotic in the fusion presence (V ^ 0, y ^ oo) 


When y ^ oo the state of the system approaches to the stationary state 
Gooiw) = limy^oo G{w; y) of eq.(^, which can be found from the equation for 
4>{w) = dGoo/dw: 

[m 2 + w(m 2 + + (’^i ~ Xw)4> = 0 . (21) 

Solutions of this equation have the form 

cl) = Ce^^{w + c)-\ ( 22 ) 

where a = X/{m 2 + v) , b = {mim 2 + miv + m 2 X)/{m 2 + , c = m 2 /{m 2 + v), 

and C is determined from the initial condition. This gives 

/ W 

(j){x)dx = 1 + C[F{w) — F{1)] , (23) 

with 

a(w+c) 

= n(,nJ-Ab ~ (-«)^”^^r(-6, -a{w + c)) , (24) 

Hi I (JJ O J 
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where F is the standard Gamma function of two arguments. 

If mi = 7712 = 0 we get 

GooH = + , (25) 

with the constant B determined by the initial condition. 

Thus, the state of the system becomes asymptotically Poisson-like: 

Mo-> 1 , Ms-> , s > 1 , 

y—*oo y—*co 

Po{y) = 1- B + Be~°- , pN -> B'^e~°- , > 1 . (26) 

y^oo I \! 

Though the moments Ms depend on the subscript s in power way the relations 
between moments reveal the correlation connected to independence of po from 
other states: 


Mi(oo) = Ba , Ms(oo) = Ba"" ^ [mi(oo)]'* . 


At B = 1 one gets the purely Poisson distribution and all correlations vanish. 

The value of the constant B which determines the asymptotical behaviour 
depends on a choice of the initial condition. Due to (^) 


B = 


1 - Go(0) 
1 - e-“ 


(27) 


Note that Go(0) = po(0) and if po(0) = Poiv) = 0) then B = {1 — e““)“^. 

If nil (or m 2 ) differs from zero the state with = 0 is the absorbing one. 
The death or annihilation processes dominate asymptotically in this case and, 
eventually, the system goes to the state without particles: 


Pn -> <5wo , Gociw) = 1 . (28) 

y^cc 

This conclusion can be done from direct analysis of the kinetic equations (0). 
after putting to zero the derivatives in the left-hand sides. 

Note that if a = X/ty ^ 1, the average number of particles in the system 
is large and its state can be described preasymptotically in thermodynamical 
approximation by eq.(^), though eventually at very large y it goes to the 
absorbing state. 


3 Account for spacial diffusion of partons 

In this section we take into account a possibility of random walk of partons 
in the transverse plane in the course of evolution. We assume that a trans¬ 
verse coordinate of a parton, b^, gets for an infinitely small interval dy an 
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(isotropic) increment dhi such that (dhi)'^ = Ddy. Note that this assumption 
on the character of the random walk of partons results in the linear form of 
the pomeron trajectory. It is a hypothesis only and other dynamical models 
for parton diffusion are conceivable 

Let us formulate now dynamical equations for the parton system. De¬ 
note the probability for the system to be in the state with N partons at 
the moment y as PNiu), and the parton distribution in the impact plane as 
Piv( 2 /; bi,..., bw). To simplify the notations we shall omit the argument y of¬ 
ten, and denote the set of N transverse coordinates as Bn = {bi,..., b^r}. All 
partons are supposed to be of the same type |^, and Pn is a symmetric function 
of coordinates normalized as 

^ J=PN{y) , (29) 

where dBN = dhi... db^, and Pn(]j) is the probability to have the A^-parton 
state at the moment y. 

The following changes of the state of the system for an interval dy are 
possible: change of position of any parton to (db^)^ = D dy] death of a parton 
with the probability mi dy; splitting of one parton into two ones with same 
coordinates with the probability A dy, and, for two partons sufficiently close 
one to another, processes of their annihilation with vanishing both partons and 
their fusion into one parton. Assuming that the parton size is small enough 
(compared to interparton distances) the annihilation and fusion probability of 
partons at points b^ and b; can be described as m 2 S{hk — b/) and uS{hk — b/) 
correspondingly Let us stress that the constants m 2 and i' are dimensional 
in contrast to the situation of sect.2. 

Thus, the evolution equations for the densities Pn have the following form: 

dpN{y]BN) ^ jjVl pN{y]BN) - {m, + X)NpN{y]BN) 
dy 

/ N^2 

dbN+iPN+i{y] Bn, b^+i) + A ^ Piv-i(y;^NV(bfc - b/) 

k,l^l 

One can consider a discrete model of parton diffusion by granulation the transverse space 
into small cells. In this case the random walk is included to the kinetic scheme on equal 
footing with splitting and fusion. The size of a cell plays the role of a cut-off for the model. 
^ If the internal quantum numbers of partons are introduced then relations between the con¬ 
stants of the RFT corresponding to given parton model, particularly between the intercept 
and the triple pomeron coupling, will be changed. 

® As it was mentioned, one can consider parton coordinates as discrete quantities in order 
to regularize divergences of the model. This results in the correspondingly smoothed <5 
functions. 
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N^2 . 

- (m 2 + X! BN)5{hk - hi) +7712 db„+ipN+2(2/; b„+i, b„+i) 

k,l=l '' 

f N 

-\-ly / dbjv+iPAr+i(j/; ^ivj bjv+i) i5(bfc bjy+i) , (30) 

fc=i 

where = EZi(^aPiv/^^L) ,a = 1,2, and Bj^^ means the set B^f, from 
which the coordinate b; is removed (Bj^^ = B^ \ bi). 

One can also define the Fourier transformed distributions 


where IC^ = {ki,... , k;^}, with normalization 

CTn( 2 /; 0 , ... , 0 ) = N\pt,{y) . 

In the momentum representation the evolution equations take the form: 
daN{y;ICr. 


(31) 

(32) 


dy 


= -D crtf{y,ICti) - {rni + X)Natf{y,ICti) 


+ mi (jN+i(y; /Cat, 0) + a E (TN-i(.y,^N\'^k + k/) 


k,l=l 


— (W 2 + u) 


E , kfc + q, k; - q) 


( 27 r) 


k,l=l 


m 2 


dq 


( 2^)2 

dq 


( 27 r) 


o-jv+2(j/;/CA,,q,-q) 

N 

^crN+i(y;/cif\kfc+ q,-q) , 


(33) 


fc=l 


where = /Cjv \ k^ and = /Cjv \ {kfc, k;} are the sets of momenta 
with the momentum kfc or, respectively, the momenta kfc and k;, removed. 

It is convenient, in the spirit of statistical mechanics, to introduce 
instead of distributions Pn(2/,^jv) a set of the multiparticle distributions 
zi,... Zs), s = 1,... ,N, which correspond to fixation of the coordinates 
of s partons and integration over coordinates of all rest: 

/f )(j/;Z,) = ^ JdB, p^iyBM^^B^lZ,) , 


(34) 
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where, as above, Zs = {zi,... , Zg}. Here the function 

Af\BN\Zs) = ^ 6(zi - bjJJ(z2 - bij ... S{zs - b^J , (35) 

where summation is taken over all sets of noncoincident indices ii, *2 ■ ■ ■ *s 
from the set {1, 2,..., A^}, is a sum of (5-function terms which fix coordinates 
of s partons from N ones. 

The total multiparticle distributions are defined by summation over states 
with different number of partons: 

My,Z,)='£fr{y:2s) . (36) 

N 

The normalization of the functions Zs) follows from eq.(p^: 

JdZsfriy,2s) = {N)sP. , (37) 

where, as above, {N)s = Nl/{N — s)! . Correspondingly, the distribution 
fs{y, Zs) is normalized to the mean factorial moment: 

JdZsfsiy Zs) =< {N)s >= Rs{y) ■ (38) 

The functions fs^^{y; Zs) are related by the reduction equation: 

nZUy2s-r) = 1^5^ / dZr+l ■ ..dzsf)r\yZs) . (39) 

Quite similarly, one can introduce the momentum multiparticle distribu¬ 
tions g)^\y; Qs), connected to the functions with Fourier transform: 

gi^^y Qs) = JdZse^^---... Zs) , 

gs{y\ Qs) = ^ 5r'(y; Qs) ■ (40) 

N 

These functions are obtained from the distributions a^iy, ICn), if to put s 
momenta equaled to qi,... ,qs, and all the rest — to zeros: 

1 

{N-s)\ 


9f\yQs) 


o-A,(gi,...,gg,0, ...,0) . 


(41) 
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The reduction relation has the form: 

Qs-r) = gf’iy; Qs-r,o,--^-,o) , 

r times 

and normalization looks as 


(42) 


ffs(j/;0, ...,0) =/is(2 /) . (43) 

The simplest way to get the kinetic equations for multiparticle distribu¬ 
tions is to reduce eq.(|^) by means of eq.(|4^) from the function to 

the function This gives 

Hj/; Q.) = -D gr(y; Q.) - (mi + X)Ngriy, Q.) 

+ mi(iV - s + l)gf+^\y, Q.) + A(iV + s - l)gf-^\y Q.) 


+ A ^ + q;) 

k,l^l 

- (m 2 + iy) f ^ 5'"^^ (j/; , q* + q, q; - q) 

k,i=i 


-2{m2 + v)J ^ 5 s+i (y; , qfc + q, -q) 

-{m 2 + i^)j J0j29i+2iy^ Qs,ci,-q)+m 2 j Q*,q,-q) 




qfc + q, -q) 



^q (iv+i) 


(271)^ 


gi+V iy^ Qs,q,-q) ■ 


(44) 


Note here, that when doing the reduction one has to consider separately the 
cases when summation variables i, k are less and greater than s, and this 
increases a number of terms with coefficients A and v. The combinatorial 
coefficients appear in the equation due to this consideration (the common factor 
{N — s)\ is omitted in the equation). 
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Part of terms cancel after summation over N (particularly terms containing 
gs+ 2 ) and one gets 

Qs) = -D gs{y] Qs) + (A - mi)sg,{y] Qs) 

s >2 , s >2 

+ A^ gs-i{y,Qi'"’^\cik + <li)-v —^ ^ 5s(?/; Qfe + q,Qi-q) 

k,l^l ^ ^ k,l^l 

-{2m2 + v)j-^^^'^gs+i{y]Qf\cik+q.,-q} , (s = l,2,...)- (45) 


The equation for fsiy^Zg) can be obtained both by Fourier transform of 
(^) or directly from the evolution equation ( ^ for using the definition (^) 
(it is necessary again to account for accurate combinatorics for cases i,k < s 
and i,k > s). As a result one gets the following equations for 

±ff\y-Zs) = DVlff\y-Z,) - (mi + A)iV/f>(y;Z,) 
yrn,{N-s + Zs) + A(7V - 1 + s)f/-^\y; Z,) 

s^2 s^2 

+ -fsZiiy, zi^^)S{zk - zi) - (m 2 + ^ /r’(2/; Zs)S{zk - zi) 

k^l — 1 k,l — l 

- 2(m2 + v)^ f''Jl\{y]Z^^\zk,Zk) - (m 2 + J^) f dufm{y;Zs,u,u) 

k=l 

+ m 2 Jdufll+^\y;Zs,u, u) + Zfc, Zfc) 

+ iyjdufll+"\y- Zs, u, u) , 




(46) 


and the set of equations for the total distribution functions /s(y; 2s)‘ 
d 


dy 


fs{y;2s)=DVifsiy;Zs)^{X-mi)sMy;Z,)^ 


s>2 


s>2 


+ h-i{y\ z2'^)5{zk - zi)-v^ fs{,y]Zs)5{zk-zi)- 


fe,Z=l 


k,l^l 


-{ 2 m 2 + i')'^fs+i{y,Z'd^\zk,Zk) , (s= 1,2,...) • 




(47) 
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Thus, we came to equations for the sets of the functions fs{y\Zs) or 
9s{y\Qs) which contain only distributions averaged over states with various 
numbers of partons. 

4 Equivalence of the parton model to reggeon field theory 


It will be shown in this section that eqs (^), are equivalent to equations 


for reggeon vertices of the reggeon field theory with a particular set of coupling 
constants — the reggeon intercept A = A—mi, the trajectory slope a' = D, the 
constant of splitting of one pomeron into two ones A, the constant of fusion 
of two pomerons into one 2 m 2 + v and the constant of pomeron scattering 
equaled to v. 

The simplest way to see this, is to make the Mellin transform of the func¬ 
tion gs{y)-. 


'OO 


Qs) , 9s{y, Qs) = ^ 


Gs{io; Qs) = / dye 


dcoe‘^yGs{io;Qs) , (48) 


T 


0 


where the integration in second formula is carried out in the complex plane w 
along the infinite contour parallel to the imaginary axis. 

One gets the following set of equations for the functions Gs{uj) 





Let the initial condition for be 


9s{y — 0; Qs) — 5ssq ■ 


(50) 


Then 9s{y', Qs) (respectively Gs(oJ] Qs) in w representation) have meaning of 
vertices of transition of sq reggeons into s reggeons. The function of free 
propagation (diffusion without splitting, annihilation and fusion) of s reggeons 
in the lo representation has the form 



(51) 
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and the renormalized propagator can be defined as 


Dsiuj; Qs) = - s(A - mi) - D ^ qlj 

If define the reduced vertices 

Qs) = [Dsito; Q,)]"' G,(w; Qs) 
the equations take the form 


(52) 


(53) 


* 

ri'^{u;-,Qs) = xY, K-A^;Qi'^^\qk + cii)Ds-i{co;Qs)- 
^ r®‘>(a;;Q('='\qfc+q,qi-q)L>s(w;Qs)- 

-{ 21712 + iy) J j^'^Tll^{uj;Qi'‘\(ik + q,-q)Ds+i{uj-,Qs) . (54) 

which corresponds to the rules of reggeon diagram techniquei and is described 
by diagrams of fig. 7. 




Figure 1: Diagrammatic representation of equations for the inclusive parton distributions 
which reproduces diagrams of RFT 


Another correspondence of eqs (^) with reggeon theory is their treatment 
as equations of “inel^ic” reggeon theooL induced by Abramovsky-Gribov- 
Kancheli cutting rulesEj (see details in ref.E^. Suppose that the contribution of 
some reggeon graph Q to the scattering amplitude is schematically represented 
as —iT^^\{—Dk}), i.e. depends on the set of reggeon propagators {D{yk, qfe)} 
integrated with required weights. Its contributicm to the total cross section is 
^tot{g) j 2 — Then (see ref. @) the inelastic discontinuity 

of this diagram, i.e. its contribution to the inelastic cross section, can be 
written as = —T^({—2Im(iZIfc)}), that is it also can be calculated 
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by means of reggeon theory but with cut propagators. Contributions of the 
diagram Q to generating functions for different inelastic characteristics also can 
be determined from this theoryE3 by changing arguments of operator T®. 

Thus, evolution of the parton distributions in rapidity is described by 
diagrams of reggeon field theory. The pomerons do not come out in this de¬ 
scription as some individual propagating objects. The evolution in rapidity 
of the s-pomeron state corresponds to the developing in time the s-particle 
inclusive distribution of partons (i.e. when momenta of s partons are fixed and 
momenta of all others are arbitrary). 

5 Interaction of colliding hadrons 

In previous sections the parton model was formulated which dynamics is equiv¬ 
alent to reggeon field theory. Parton distributions at different y given by this 
model can be associated with the stationary structure of parton distributions 
in rapidity for a fast moving hadron (i.e. with its Fock wave function). 

We define in this section the “interaction operator” of two parton systems, 
which allows to calculate the amplitude for interactions of two fast colliding 
hadrons. The form of this operator has to meet the requirement of correspon¬ 
dence with calculations according to the RFT rules. 

Let two hadrons h and h move one toward another with rapidities y and y 
in some Lorentz frame at impact parameter b. The essential feature required of 
the interaction operator is that Lorentz invariance of the scattering amplitude, 
that is, its independence from a choice of particular Lorentz frame, should be 
maintained. In other words, the amplitude should be a function of the total 
rapidity Y = y + y only ^ , but not of the individual values of y and y. 

The hadron states are specified by sets of parton distributions {fs{y;Zs)} 
and {fs{y\Zs + b)}, where designation Zs + h means a shift of all coordi¬ 
nates Zq to the vector b in the impact plane. Denote for brevity these sets as 
vectors in space of states |iF(j/;0)) and |iF(y;b)). In the momentum represen¬ 
tation we shall use notation \G{y', 0)) and \Giy\ b)) for the sets {gs{y', Qs)} and 

{9s{y; Qs)exp(ibX;i Qa)}- 

Let the total cross section is determined as an integral of function T{y, y\h) 
over all impact parameter values: 



(55) 


^ Though the hadrons move in the opposite directions, we define both y and y as positive. 
So the total rapidity which is the rapidity value of one hadron in the rest frame of other is 
written as a sum Y = y y. 



Probabilistic Model of Reggeon Field Theory 19 


We introduce the “interaction operator” T by the following equation: 
T{y,y;b) = {:F{y;h)\f\T{y;0)) = 

= ^ - J dZsdZsfs{y-,Zs)fs{y,Zs)S^^H2s - i. - b) , (56) 


or, in momentum space, 


T{y,y-,b) = {g{y-h)\T\g{y;0)) = 


= E 


(- 1 )^ 


5 ! 


dQs 

(27r)2 


jb 1 qc 


gs{y-,Qs)9s{y,Qs) ■ (57) 


The parameter e plays a role of the squared size of a parton. As we mentioned 
in sect.3 in order to regularize some singular expressions the function 
should be smeared out over this range. 

Let us show that the definition (^), ( p7|) of the operator T combined with 
particular relation between constants provides the dependence of the function 
T(j/, y; b) only from the single variable Y = y + y. In the momentum represen¬ 
tation the evolution equations (®) can be symbolically written as 


-\g{y-h)) = H\g{y;h)) , 
dy 

where the form of the operator H follows from eq. (|4^) . Its solution is 


(58) 


|5(y;b)) = exp 


Hy 


1 ^(0; b)) . 


(59) 


i.e. the evolution operator exp 
erator. Therefore 


Hy has a meaning of the Lorentz boost op- 


T{y,m = {g(j)-Me^"yfe^y\G{S)M 


(60) 


{H* is the transpose of H). If the relation 


(x|i7‘T|</>) = {x\fH\4>) (61) 

holds for arbitrary states \4> >, \x >, then the required property will be fulfilled: 

(X|e^‘"fe*"|<(.) = (x|fe"(^+^)|(^) . (62) 
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Let us prove the relation (|^. The first, second, and fourth terms of the 
matrix H (see eq.(^5[)) are diagonal in the parton numbers s, and the relation 
( 1^ ) is evident. The third and fifth terms (denote them as H 12 and H 21 ) 
correspond to splitting and fusion of reggeons and their matrix elements have 
to be conjugated to each other, i.e. {x\TH 2 i\(l)) = {x\Hl 2 T\ 4 >) and vice versa. 
Indeed, the matrix element of TH 21 is written as 


{x\TH2i\(I)) = - {27712 + 


\s—1 


(_l)s-.e^ 


dQs dq 


S =1 


(27r)^® (27r)^ 


XsiQs) 


^ / 1 ^ ® ^ 
^'p8+i{Qi^\qk + q,-q) = {2m2 + u)^——^ 




S =1 


(s-1)! 


X J j^Xs{Qs)(ps+i{Qi'"\qk + q,-q) , (63) 


with regard that the sum over k contains s identical terms. In the matrix 
element of TH 12 summation over parton number starts with s = 2. Hence, 
after substitution s = s + 1 and taking into account that the sum over /c, I 
contains s(s + 1) identical terms one gets 


/_ 1 'js —1 

ixiHum ={7j7\TH,2\x) = 

S=2 

S 

X Y xs-i{Qf''\<ik + qi) 

k,l^l 

f dQs^i _ 


dQs 

(27r)2S 


(t>s{Qs) 




(-l)^e« 

(s-1)! 


+ Cls+l)?!>s+l (Qs+l) • 


(64) 


After change of integration variables the both matrix elements are the same 
provided 


2 m2 + z/ = A e . (65) 

In the case of the 0-dimensional model constants A, m 2 and 1 / are dimensionless, 
the parameter e in the definition of the interaction operator is not needed and 
the requirement of Lorentz invariance gives the relation 

2m2 + 17 = X ■ (66) 

The conditions (^5|), (|^) look quite natural from the point of view of RFT 
(equality of three-pomeron vertices). 
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Let us discuss briefly relation between diagrams representing a history 
of development of parton cascades and the RFT diagrams. Fig.2a gives an 
example of evolution of parton distributions for the colliding hadrons with 
rapidity. In this example two partons from each hadron are linked in some 
intermediate Lorentz frame (s = 2 in eqs (^), (|^). The cascade branches 
with “dead ends”, i.e. terminated in noninteracting partons (thin lines in 
fig.2a), result in renormalization of distribution functions (see sect.4) after 
summation over all possibilities of the “dead-end” branches. The renormalized 
propagators are showed by thick lines at the figure and form a skeleton of the 
parton network. Each skeleton corresponds to a particular order of reggeon 
perturbation theory in number of splittings and fusions of renormalized reggeon 
propagators. Different ways of linkage of partons from interacting hadrons give 
different skeletons, i.e. different reggeon diagrams. Skeleton can be imaginated 
pictorially as a set of lines in the cascade network with non-zero electric current 
flow if to ascribe definite potentials to the sets of initial points of the projectile 
and target. 

Note that the interaction operator ( ^^ doesn’t have a sense of probabil¬ 
ity of interaction for two parton systems because it takes into account only 
single interactions of partons. Similar situation one can see in desertion of 
nucleus-nucleus scattering in so called Czyz-Maximon approximation til which 
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includes only first-order interactions of every nucleon. It describes correctly 
the scattering amplitude only for low nuclear densities, and account for all 
orders of nucleon interactions is necessary to get unitary amplitude for high 
densities. However, in case of parton model, introduction of higher order terms 
into interaction operator would violate Lorentz invariance of the amplitude. 


6 Examples 

6.1 Eikonal approximation 

Consider first hadron scattering in the model including parton diffusion but 
without any parton transformations (mi = m 2 = A = = 0). In this approxi¬ 

mation each parton walks randomly and independently in the transverse plane, 
therefore, if there were no correlations in the initial parton distributions for 
both hadrons there will be no correlations in the course of further evolution: 

S S 

fs{y; 2s) = fi(y; Za) , fs(y] 2s) = fi(y; Za) ■ (67) 

a—1 a—1 

The amplitude of hadron interaction at the impact parameter b is determined 
by eq.(^) and coincides with the eikonal approximation formula (see fig.3a) 

Tie^k) ^ ^ x(r, 6) = 1 - , (68) 

where Y = y + y, and so called eikonal function is equal to 

x(Y,b) = e J dz/i(y,z)/i(y,z - b) = e j j^^e^'^^g(y,q)g{y,q) . (69) 


If the initial distributions in hadrons were of the d-function type, then 
one-particle distribution at finite y is determined by random walk only: 

/i(y, 2 ) = , gi(y,q) = (70) 


and similarly for the hadron h. This gives 


X(Y, b) 


e 

AttYD^ 


4YD 


(71) 


If the initial distribution is given by some function r(z) with Fourier transform 


s(q) = 


dzE^^r(z) , 


(72) 
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interpretation of the eikonal approximation for supercritical pomeron 

then the eikonal is determined by formula 

X{Y, 6) = e y ^ ■ (73) 

It is seen that eq.(^ corresponds to exchange of the regge pole with intercept 
equaled to 1, trajectory slope a' — D and residue s{q)s{q). 

Often the formulae of the eikonal approximation are applied in the case of 
supercritical pomeron with intercept larger than 1. It seems to be inconsistent 
in the model under discussion because in the parton picture the supercritical 
regime means the growth of parton number when rapidity changes and this can 
occur only due to parton splitting in the course of evolution, i.e. at nonzero 
value of the constant A. The eikonal approximation in the presence of split¬ 
ting means that among soft partons originated due to evolution of each initial 
partons only one may interact with a target (see fig.3b). The more consistent 
consideration of parton splitting allowing interaction with the target for any 
parton is carried out in the next example. However, the analytic formulae can 
be obtained only if parton diffusion is disregarded. 

6.2 The Q-dimensional Schwimmer model 

The Schwimmer modelB was developed for description of inelastic interaction of 
fast hadron with heavy nucleus and corresponds to summing enhanced reggeon 
diagrams in the tree approximation (see fig.4). If a radius of the nucleus is 
large enough (i?^ ^ ofY), then parton diffusion can be disregarded and the 
problem is equivalent to the 0-dimensional model of sect.2a (we assume that 
A 7 ^ 0, mi 7 ^ 0, m 2 = v = Q). 

Choose the initial state of the projectile hadron as containing one parton with 
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Figure 4: Reggeon diagrams corresponding to the Schwimmer approximation 


a probability gp 

= 9p , p'o'’ = ^- 9p , P^N at iV > 2 , 

G^P\w,Q) = l + gp{w-l) , (74) 

and the initial parton distribution for the target in the Poisson form 

N 

P%\Q) = %e-P\ , gW{0) = {g,y. (75) 

The fusion is absent in this model, = 0, so the Lorentz-invariant consideration 
is impossible. The interaction amplitude calculated in the laboratory frame 
has according to the 0-dimensional version of the formula (^) the following 
form: 


/_-I '\s—1 

TschAY) = E M— 


Si 


= 1 - - gt;Y) = 


9p9t^ 


AY 


l+5^(gAV_i) 


(76) 


We used the definition (j^) for the generating function G^p^ and its explicit form 
(0), disl) as applied to the initial condition (^: G^p\w; Y) = gpW{Y), with 
W{Y) from eq. 
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When one estimates the cross section of hadron interaction with a nucleus 
A one has to choose gt to be proportional to the nuclear density p^(b): 

gt=gNPA{h). (77) 

For heavy nucleus (A ^ 1) in the approximation of constant density, 

Pa = -^0{Ra - |b|) , Ra = , (78) 

irRi 


we come to 


T (Y h^ - 3 p 9 nA^/^ exp(Ay) 6»(i?^ - |b|) 

, Aa.Ai/Vexp(Ay)-l 


(79) 


This formula coincides with the original Schwimmer formula for amplitude of 
hadron-nucleus scattering in the approximation of very large nucleus. When 
Y ^ oo the amplitude goes to the constant limit gpA/X = gp{l — mi/X) < 1. 
Note that in general reggeon field theory the constants gp, A and A may be 
arbitrary. 

It is worth to stress that though the formula ( [7^ coincides in form with the 
Schwimmer formula obtained in RFT, the meaning of the constants entering 
these formulae is different, and even their dimensions differ: in eq.(f^) the 
constants gp and A are dimensionless, and g^i has dimension of length squared, 
while in reggeon field theory all these constants have dimension of length. 

It is easy to formulate a generalization of the Schwimmer approximation 
in this approach. If the initial parton distribution in the projectile has the 
Poisson form, i.e. 


= (PpY , G^P\w,0) = ex.p[gpw] , G^pYw.Y) = ex.p[gpw{Y)] , (80) 


one gets the “eikonalized” Schwimmer formula 

(^) = 1 - exp[T,,i,^ (F)] . (81) 


Note that analysis of interaction of two heavy nuclei requires even in the 
tree approximation taking into consideration not only parton splitting but their 
fusion also. Moreover, as it was emphasized already, the Lorentz-invariant 
parton consideration is possible only in presence both splitting and fusion of 
partons provided the corresponding constants are related. 
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6.3 Asymptotics of the 0-dimensional model with account for parton fusion 

We consider in this section the 0-dimensional model including splitting and 
fusion of partons for the case mi = m 2 = 0. Let us demonstrate the Lorentz 
invariance of the parton picture for asymptotically large rapidity value Y, when 
a stationary regime sets in and the interaction amplitude does not depend on 
Y at all. Indeed, in this case for any frame choice one of hadron rapidities (say, 
y) is large. Hence the projectile hadron is in asymptotical stationary state (see 
sect. 2.b) 


= a = Xfn = l, (82) 

1 — e ^ 

and the state of a target with rapidity y is given by the generating function 
G(*)(y, w). It is essential that the value of this function at w = 0 does not 
depend on y and is equal to Gq^(O). 

Therefore for a = 1 the interaction amplitude does not depend on rapidity 
in any Lorentz frame and equals 


^ / _ 1 1 

TUV) = E |l - G«(y, 1 - a) 


Si 


[1_gW(0)][1-G(,'^(0)] 


(i)/ 


1 — e 


-1 


(83) 


5. Two-dimensional diffusion 

In presence of parton diffusion it is possible to distinguish in the parton struc¬ 
ture of the fast hadron two characteristic regions. At the periphery where 
parton density is small the parton fusion is insignificant and only diffusion and 
splitting processes occur. Parton distribution in this region behaves like the 
Green function of diffusion equation complemented with exponential increase 
due to splitting: 


p{periph) 

Jl 




= gAy 

Airy 


(84) 


In the central region the parton fusion is essential and the density comes to 
the stationary value 




( 85 ) 
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Transition from one regime to another corresponds to equality of the central 
and peripheral density values, i.e. to distances of order 

R{y) « 2VXDy . (86) 

At larger distances the density decreases exponentially. Thus, the fast hadron 
represents to be a disc of constant density of radius increasing proportionally 
to its rapidity: 


fi{y,z):^^9[R{y)-z] . (87) 

Obviously the interaction cross section of two such discs corresponds to the 
Froissart regime: a{Y) ^ const • Y^, where Y = y + y. In order to calculate 
the constant value, i.e. the degree of blackness, it is necessary to account 
for contributions of all multi-particle densities fs{y,Zs), fs{y,Zs). It is not 
the aim of this paper to carry out complete analysis of this problem. Let us 
consider as an illustration only the simplest case when the initial state of one 
of hadrons (say, h) contains only one parton, that is fi{y = 0,z) = 6{b — z), 
fs{y = 0, Zg) = 0, s > 2. Let us use the Lorentz invariance of the problem and 
apply the formula ( |^ ) in the rest frame of the hadron h. Then in the formula 
only the first term remains: 

T{y,y;b) = e J dzfi{Y; z)f i{0; z)6{z - z - b) 

^ — 9{R{Y) - b) = 9{R{Y) - b) , (88) 

u 

which corresponds due to relation (^) to scattering on the absolutely black 
disk of radius R{Y). 

For more complicated cases (many-parton initial states) the situations 
turnes out to be the same — the asymtotics of the amplitude corresponds 
to the scattering off a black disc. In ref. Ilj Kancheli discussed the problems 
related to the Lorentz invariance of interaction of parton systems for the case 
of the gray disc (with finite transparency). In present model this problem does 
not appear, and both Lorentz invariance and blackness of the amplitude are 
provided by relation (^). 

7 Unitarity 

One of the key problems of RFT is s-channel unitarity of the scattering ampli¬ 
tude. As can be seen, the unitarity of the amplitude calculated for the model 
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under consideration depends not only from the parton dynamics but also from 
the initial conditions, i.e. from the form of the parton sources for each hadron. 

Let us discuss under which conditions the amplitude T{Y, b) can be inter¬ 
preted as a probability of interaction at fixed impact parameter that is satisfies 
to the relation 


0<T{Y,b)<l. (89) 

Consider first the case with d = 0 and note that if both sets of moments 
fj-siy), Bs{Y - y) in equation 


°° 1 

TiY,b)^Y^ -j- Hs{y)fLs{Y - y) (90) 

s=i 

would be of Poisson form, 

= a'* , = d'* , (91) 

the amplitude has a simple form satisfying to relations (^^: 

(F, b) = l- exp(ad) . (92) 

This observation makes useful so called Poisson representationlli for prob¬ 
abilities. It represents (for the simplest case) probabilities pwiv) as superpo¬ 
sitions of Poisson distributions Pat (a): 


PN{y) = J daf{a,y) j^e “| = (Pv(a)) , 

where we used notation (...) for averaging with weight /(a, y). 
The generating function G{w; y) can be written as 

G{w] y) = J da f{a, y) exp[(u; - l)a] . 


(93) 


(94) 


The master equation (^ gives for f{a,y) the equation of Fokker-Planck type 


df{a,y) 

dy 




(95) 


where functions 


A{a) = (A — mi)a — + 2m2)a^ , B{a) = 2Xa — 2{v + m 2 )a^ , 


(96) 
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describe driving term and diffusion coefficient correspondingly. The diffusion 
coefficient B is non-negative for 0 < a < A/(z^ -I- m 2 ). The condition (||) 
provides normalization of f{a,y): 

J daf{a,y) = 1 . (97) 

Equation is equivalent 0 to the differential stochastic equation 

da = A{a)dy -sjB{a)dW(y) , (98) 

where dW(y) is a stochastic differential corresponding to the Wiener process. 

Suppose for simplicity that mi = m 2 = 0. One can see from (^) that if 
evolution starts at ?/ = 0 at a value within interval (0, o), a = A/zz, then the 
first term will drive the system to the right end of the interval a = X/v). At 
this point both the driving term and diffusion coefficient vanish and, therefore, 
it is a stationary point. Any initial a-distribution will be concentrated at this 
point when y —> 00 . Thus the asymptotic distribution is Poisson distribution 
PAr(a), in correspondence with analysis of section 2. 

Because of positivity of the diffusion coefficient B{a) at the interval (0, a), 
the distribution /(a, y) is positive and normalizable in the course of evolution, 
i.e. it can be considered as a probability distribution. If both f(a,y) and 
/(a, y) are of class of probability distributions, then the interaction amplitude 
T{Y, b) is given by averaging of equation (^2|) over a and d, 

T{Y,b) = {{T^^’^\Y,b))) = J daf{a,y) J dd/(d, y) [1 - exp(ad)] , (99) 

and evidently satisfies the unitarity condition (|8^). Consideration of two- 
dimensional case can be carried out similarly. 

8 Conclusions 

The parton stochastic model discussed in this paper is in exact corre¬ 
spondence with reggeon field theory. Only special reggeon field theories with 
particular set of pomeron vertices allow parton interpretation of this sort. 
The model reproduces results of main approximations for RFT. In the the¬ 
ory with zero number of transverse dimensions, d = 0, it is possible to analyze 
quantitatively an asymptotical regime. It is equivalent to summing all dia¬ 
grams of 0-dimensional RFT with all loops. For theory with d = 2 qualitative 
and numerical consideration is possible. Because of its classical and stochas¬ 
tic nature the model allows numerical simulation by means of simple Monte 
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Carlo algorithm. Direct analogy with kinetic theory of chemical reactions gives 
good prospects for application of thermodynamical methods and of theory of 
stochastic equations. 

The model allows a number of generalizations, in particular to parton 
systems with internal quantum numbers and of various spacial scales. 
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